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Abstract 

We study the full counting statistics (FCS) in a single-molecule magnet 
(SMM) with finite Coulomb interaction U. For finite U the FCS, differing 
from U — > oo, shows a symmetric gate- voltage-dependence when the coupling 
strengths with two electrodes are interchanged, which can be observed experi- 
mentally just by reversing the bias-voltage. Moreover, we find that the effect of 
finite U on shot noise depends on the internal level structure of the SMM and 
the coupling asymmetry of the SMM with two electrodes as well. When the 
coupling of the SMM with the incident-electrode is stronger than that with the 
outgoing-electrode, the super-Poissonian shot noise in the sequential tunneling 
regime appears under relatively small gate- voltage and relatively large finite U , 
and dose not for U — > oo; while it occurs at relatively large gate- voltage for 
the opposite coupling case. The formation mechanism of super-Poissonian shot 
noise can be qualitatively attributed to the competition between fast and slow 
transport channels. 
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1. Introduction 



Electronic transport through a single- molecule magnet (SMM) has been in 



tensively studied both experimentally [1] and theoretically [! 
stimulated by the prospect of a new generation of molecule-based electronic and 
spintronic devices [10j. Recently, the current fluctuation in electron transport 
through single-molecule magnet or molecular junction has been attracting much 
interest^ la [ill HiQjQjQl owing to its allowing one to identify the internal 
level structure of the transport system 16|| and to access information of electron 
correlation that can not be contained in the differential conductance and the 
average current [l3|- These studies were mainly focused on the infinite Coulomb 
interaction regime In fact, Coulomb interaction is usually finite in a re- 

alistic mesoscopic system and thus one should consider the effects of the finite 
Coulomb interaction on the current correlation. In most mesoscopic systems, 
the negative correlation induced by Coulomb interaction may impose a time 
delay between two consecutive electron transfers and lead to the suppression of 
shot noise so that sub-Poissonian shot noise occurs. For example, in symmet- 
rical double-barrier i unctions [lq and in nondegenerate diffusive conductors [l^] 
the 1/2 and 1/3 suppression Fano- factors have been found, respectively. How- 
ever, in the presence of a strong nonlinearity of the I-V characteristics, the 
Coulomb interaction can yield a positive correlation and enhance the noise even 



to be super-Poissonian[20, 21]. This phenomenon was first discovered in double- 



barrier tunneling diodes in the negative differential conductance (NDC) regime 
and a Fano factor up to 6.6 was observed [20J]. In the NDC regime, Coulomb 
interaction and the density-shape of states in the well introduce positive correla- 
tion between consecutive current pulses, which leads to a super-Poissonian shot 
noise. In general, the effect of Coulomb interaction on the shot noise is more 
complicated. For instance, in mesoscopic coherent conductors (at sufficiently 
large voltages) Coulomb interaction decreases the shot noise at low transmis- 
sions and increases it at high transmissions [22J. Shot- noise measurements in 
mesoscopic devices can also provide information about the effective charge e* of 
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the current-carrying particles. For a quantum-dot system in the Kondo regime, 
the simultaneous presence of one- and two-quasiparticle scattering results in a 
universal average charge e* = 5/3e[23j. However, as shown in Ref. [24j], the 
Coulomb interaction remarkably influences the effective backscattering charge 
of current-carrying particles via a correction factor (e* = 5/3eF) which is less 
than unity. Furthermore, the effect of Coulomb interaction on the shot noise can 
be employed to reveal important information of the energy profile of nonequilib- 
rium carriers injected from an emitter contact, but which can not be obtained 
from shot-noise measurements in the absence of Coulomb interactions [2g. I n 
the present SMM system with finite Coulomb energy U, an electron from one 
lead tunnels into the SMM and then leaves the SMM through the other lead via 
two kinds of transition processes: (i) between the singly-occupied and empty 
states, (ii) between doubly-occupied and the singly-occupied states which does 
not occur under the infinite Coulomb interaction. Therefore, it is significant to 
study the effect of finite Coulomb interaction on the current fluctuation in the 
SMM system. 

An alternative way to investigate the current fluctuation, known as the full 
counting statistics (FCS), has been proposed in the pioneering work by Levitov 



et al. 



26|| . The method yields not only shot- noise power but also all the statis- 



tical cumulants of the number of transferred charges. The transport through 
mesoscopic devices is fully described by the FCS, which may provide the full in- 
formation about the probability distribution P (n, t) of transferring n electrons 
between electrode and SMM during a time interval t. The FCS is obtained 
from the cumulant generating function (CGF) that is related to the probability 
distribution byQ 

e -F(x) =^p( n ,i)e mx , (f) 

n 

where x is the counting field. All cumulants of the current can be obtained 
from the CGF by performing derivatives with respect to the counting field 

C\ = — (— id x ) k F (x) . In the long-time limit, the first three cumulants 

x=o 

are directly related to the transport characteristics. For example, the first- 
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order cumulant (the peak position of the distribution of transferred-electron 
number) C\ — n gives the average current (J) = eC\/t. The zero- frequency 
shot noise is related to the second-order cumulant (the peak-width of the dis- 
tribution) S = 2e 2 C 2 /t = 2e 2 (v? - n 2 ^ /t. The third cumulant C 3 = (n - nf 
characterizes the skewness of the distribution. Here, (• • • ) = J2n (' " ) P ( n > 
In general, the shot noise and the skewness are commonly represented by the 
Fano factor F2 = C2/C1 and F3 = C3/C1, respectively. 

Since the electron-electron interaction may bring correlations and entangle- 
ment of electron states, the shot noise in the mesoscopic system with Coulomb 



interaction has attracted the significant attention [22|, |24|. The study of FCS of 
interacting electrons in mesoscopic systems has become a challenging subject 
of great interest [27[. In this letter, we study the FCS of electron transport 
through SMM weakly coupled to two metallic electrodes, here the first three 
cumulants are given in terms of numerical calculation. In the model considered 
here two electrodes are regarded as non-interacting Fermi gases, while the cen- 
tral SMM is treated as a multi-level system with finite Coulomb interaction. We 
found that the effect of finite Coulomb interaction on the shot noise depends 
not only on the internal level structure of the SMM, but also on the left-right 
asymmetry of the SMM-electrode coupling. In particular, our analytical re- 
sults indicate that for finite Coulomb interaction the FCS, which is different 
from the case of U 00, shows a symmetrical gate- voltage-dependence when 
both the intensities of the SMM coupling to the left and right electrodes are 
interchanged, which originates from the both symmetries of the effective chan- 
nel energy levels and the probability distribution. Moreover, we also observed 
super-Poissonian noise for symmetric coupling situation. The paper is organized 
as follows. In Sec. II, we introduce the SMM system and outline the procedure 
to obtain the FCS formalism based on a particle-number-resolved quantum mas- 
ter equation. The numerical results are discussed in Sec. Ill, where we analyze 
the occurrence-mechanism of super-Poissonian noise and discuss the effects of 
Coulomb interaction, the left-right asymmetry of SMM-electrode coupling, and 
the applied gate voltage on the super-Poissonian noise. Finally, in Sec. IV we 
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summarize the work. 

2. MODEL AND FORMALISM 

A magnetic molecule coupled to two metallic electrodes L (left) and R (right) 
is described by the Hamiltonian[5| H to tai = H mo i+HL ea d s +HT- We assume that 
the SMM-electrode coupling is sufficiently weak so that the electron transport is 
dominated by sequential tunneling through a single molecular level with on-site 
energy Sd- The molecular Hamiltonian is given by 



Here the first two terms depict the lowest unoccupied molecular orbital (LUMO), 
n = d\.d^ + djrfj. is the number operator, where S a (d a ) creates (annihilates) 
an electron with spin a and energy Sd (which can be tuned by applying a gate 
voltage V g ) in the LUMO. U is the Coulomb interaction between two electrons 
in the LUMO. The third term describes the exchange coupling between electron 
spin s in the LUMO and the giant spin S, the electronic spin operator s = 
^2 aa i d\ (0*0-0-') d<r', where a = (a x , <r y , <r z ) denotes the vector of Pauli matrices. 
The forth term is the anisotropy energy of the magnetic molecule whose easy- 
axis is Z-eods (K > 0). The last term denotes Zeeman splitting. For simplicity, 
we assume an external magnetic field B is applied along the easy axis of the 



The relaxation in the electrodes is assumed to be sufficiently fast so that their 
electron distributions can be described by equilibrium Fermi functions. The 
electrodes are modeled as non-interacting Fermi gases and the corresponding 
Hamiltonian 



where a akcr (a Q k<x) creates (annihilates) an electron with energy e Q ko-, momen- 
tum k and spin a in a (a = L, R) electrode. The electron tunneling between 
the LUMO and the electrodes is described by 



H mol = {e d - eVg)n + -h{h - 1) - J s ■ S - KS 2 Z - B{s z + S z ). (2) 



SMM. 




(3) 




(4) 
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According to Ref. Q the eigenstates of an isolated SMM have four branches 
and may be denoted as \n, m) v where n (n = 0, 1, 2) is the electron number in 
the molecule orbital, and m is the quantum number for the Z-component of 
the total spin. The index v{= ±) appears only when n = 1. In term of the 
electron state |i, j) LUM q — 0, y, J,) in molecular orbital and the local spin 
state \m) GS (to G [—5,5]) the empty-branch and doubly-occupied states may 
be expressed as 

|0, m) = |0, 0) LUMO \m) GS ,me[-S,S], (5) 

and 

|2, to) = |t, i) LUMO \m} GS , m E [-5, 5] , (6) 
respectively, and the two singly-occupied branches are |1, ± (5 + 1/2)) = \a±) ±5) 



(a + =t,a_ =1) for m = ±(S + 1/2), and 



|l,m) =a±\\) LUMO 

with 



m -\) + b m\i) L UM<> "' ., / 
z ' GS z / 



, (7) 



GS 



at 



J^5 (5+1) -to 2 + 1/4 
2^/AE (m)y/2AE (to) + (2X - J) to' 
, ^/2AE (m) + (2K - J)m 

m ~ T 2y/AE (to) ' 

for -5+1/2 < m < 5—1/2, where AE (to) = K (K - J) to 2 + (J/4) 2 (25+1) 
The corresponding energy eigenvalues of molecular eigenstates are given by[5[ 

e (0, m) = -Km 2 - Bm, (8) 

e(2,TO) = 2(£ d - e y g ) + C/-XTO 2 -BTO, (9) 



e* (1, m) = £d - eV^ - 5m + - - K \ mT + -j±AE (to) . (10) 

Here, for |1, ± (5 + 1/2)} the upper (lower) sign applies if K — J/2 is positive 
(negative). It is found from Eqs. (15|)~ ([T(H) that the energy eigenvalues of |0,m) 
is independent of gate voltage while those of |l,m) ± and |2,to) depend on 
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gate voltage, so that the internal level structure can be tuned by an applied 
gate voltage. For the finite Coulomb interaction case one need consider the 
contribution of the doubly-occupied states |2,m) to the transport. 

In sequential tunneling regime the transitions between the SMM and the 
electrodes are well described by quantum master equation of a reduced density 
matrix spanned by the eigenstates of the SMM. The FCS of electron transport 
through the SMM may be implemented with the help of the particle-number- 
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resolved master equation for the reduced density matrix which is given by 

(t) = -iCpW (t) - \^lp (n) (t) (11) 

under the second order Born approximation and Markovian approximation, with 
Liouvillian superoperator £(•••) = [H mo i, (•••)] and 

KpW (t) = [44~V n) (*) + p (n) (*) 4 +) 4 - 4"> (n) (*) 4 

-4 P (") (*) 4+> - 4:> (n ^ 1) (*) 4 - d > [n+1) (*) a rI\ + H ^ 12 ) 

Here ^ — 5Z a =£,ii^ a M , ^4q^ = r a n^ (— £) d^, = f a , n a = 1 — f a (f a 
is the Fermi function of the electrode a), T a — 2irg a {taf 1 and g a (a = R,L) 
denotes the density of states of the metallic electrodes. p( n > (t) is the reduced 
density matrix of the SMM conditioned by the electron numbers arriving at the 
right electrode up to time t. Throughout this work, we set e = K = 1. The 
CGF connects with the particle-number-resolvcd density matrix by defining 
S(x,t) = E„P (n) {t)e inx . Evidently, we have e~ F ^ =Tr[S (x,t)\, where the 
trace is evaluated over the eigenstates of the SMM. Since Eq. (jlll) has the 
following form 

p(n) = Ap {n) + Cp (n+1) + £)p (n-l) ) ^3) 

then S (%, t) satisfies 

S = AS + e~ lx CS + e ix DS = £ x S, (14) 

where S is a column matrix, and A, C and D are three square matrices. Here, 
the procedure for calculating the specific form of C x is given in detail in the 
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Appendix. In the low frequency limit, the counting time (i.e., the time of 
measurement) is much lon ger than the time of tunneling through the SMM. In 
this case, F (x) is given by 
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31] 



F( X ) = -\i (xK 



(15) 



where Ai (x) is the eigenvalue of C x which goes to zero for x ~ ► 0. According 
to the definition of the cumulants one can express Ax (x) as 



A 



OO 



fe=l 



fc! 



(16) 



Inserting Eq. (|16p into \C X — \\ (%) I\ — and expanding this determinant in 
series of (ix) k , one can calculate Ck by setting the coefficients of i\ equal to 
zero. 



3. NUMERICAL RESULTS AND DISCUSSION 

In the following we study the effect of finite Coulomb interaction on the 
FCS of electronic transport through the SMM weakly coupled to two metallic 
electrodes. We assume the bias voltage (VJ, — (J-l — A*r) is symmetrically entirely 
dropped at the SMM-electrode tunnel junctions, which implies that the levels of 
the SMM are independent of the applied bias voltage even if the couplings are 
not symmetric. The parameters of the SMM are chosen asQ S = 2, Sd = 200r, 
J = 100I\ K = 40I\ B = 80r and k B T = 10T, where V is the typical rate for 
the tunneling of electrons between the SMM and electrode. Here, the validity of 
the FCS formalism in this work deserves some discussions. In the SMM system, 
the relaxation rate of transition from a given molecular eigenstate \n, m)of the 
SMM to the neighboring eigenstates \n, m ± 1) of the same multiplet can be 
written as[8j 

■p\n,m)— v|n,m±l) ^ 1 /-,-\ 

SMM ~~ T " \ 77 rp ' 

TSMM 1 -f e -(,e|n,m)-6|»,m±l)J/«BT 

where tsmm denotes the SMM's spin-relaxation time, and €\ ntTn \ the eigenvalue 
of the molecular state \n, m). For a typical SMM at low temperature (about the 
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order of IK ), e.g., Fes, TSMM is of the order of 10 -6 s[32[, whereas the tunneling 
time of electron through the SMM To of the order of 10 -9 s\^. This means that 
an electron, which from incident-electrode tunnels into the SMM, has enough 
time to tunnel out the SMM before the corresponding molecular eigenstate re- 
laxes to its neighboring eigenstates of the same multiplet since tsmm ^ T o- 
Therefore, the effect of the SMM's spin-relaxation processes on the FCS can be 
neglected. The thermal energy corresponding to the temperature considered is 
much smaller than energy barrier so that thermally activated transitions above 
the barrier can be neglected. Furthermore, we assume the external magnetic 
field is applied along the SMM's easy axis and transverse anisotropy is small 
enough relative to the easy-axis anisotropy so that the magnetic quantum tun- 
neling can also be neglected. 

In the present work, we only study the transport above the sequential tun- 
neling threshold, i.e., VJ, > 2e se , where e se is the energy difference between the 
ground state with charge N and the first excited states N — 1 34] . In this regime, 



the inelastic sequential tunneling process is dominant, thus electrons have suffi- 
cient energy to overcome the Coulomb blockade and tunnel sequentially through 
the SMM. A special emphasis is the effects of the strongly asymmetric coupling 
to two electrodes and gate voltage on super-Poissonian noise. The sequen- 
tial tunneling requires a change of the electron number by An = ±1 and the 
magnetic quantum number by Am = ±1/2, so that the sequential tunneling 
threshold depends on the effective channel energy levels in bias voltage window 
which are given by 



ef (f) = e ± (l, m + i) - e(0, m) = e d - V g + E ± (m + 1), ( I s ) 



4(1) = e(2,m)-e ± (l,m+i) = e d -V g + U - E ± (m+^), (1-)) 
ef (|) = e(2, m) - e ± (l, m - \) = e d - V g + U - E±(m - \), (20) 



2 v 2' 

4 {{) = e ± (l, m - i) - e(0, m) = e d - V g + E ± (m - ~), (21.) 
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where E ± (m + l/2) = - (m + 1/2) K - B/2 + J/4 ± AE(m + 1/2) and E ± (m~ 
1/2) = (to- l/2)K+B/2+J/4±AE(m-l/2). The electron with spin <j(=t,i) 
can be transferred by the effective channel energy levels ef(a)(i = 1, 2, 3, 4). The 
presence of finite Coulomb interaction may induce that the effective channel 
energy levels move symmetrically in the bias voltage window with increasing 
the bias voltage. From Eqs. (fT51) - (j2"Tj) . one may find |ef(t)| = \ef (4-)| and 
| e 3 (t)| = | e 4 (4-) I when V g = e d + U/2, which indicates the effective channel 
energy levels ef (J.) and (t), ef (4-) and (f) enter synchronously the bias 
voltage window with increasing the bias voltage. When the gate voltage has a 
departure from V g = e d + U/2, for example, V g = e d + U/2 ± AV, the effective 
channel energy levels still have similar symmetry, 

\4(t)\ (V g = e d + U/2 ± AV) = \4 (|)| {V g - e d + U/2 T AV), (22) 

\4U)\ {V g = (d + U/2 ± AV) = \ef (t)| (V g - e d + U/2 =p AV). (23) 

For the case of U = 100r, the current as a function of the bias voltage for 
r.L/r_R = 10 and 0.1 is shown in Figs. 1(a) and (d). Both currents for the gate 
voltage V g = 100r and 400r (V g = 200r and 300r) have same sequential tun- 
neling threshold, which corresponds to AV = 150r (AV = 50r). In particular, 
we find that the finite Coulomb interaction may induce that the FCS shows the 
same bias- voltage-dependence under different external conditions. Figures 1(a) 
and (d), (b) and (e), and (c) and (f) show the average current, shot noise and 
skewness for the two different strong asymmetric couplings (Tl/Tr = 10 and 
0.1) in the presence of finite Coulomb interaction (U — 100r), respectively. It is 
interesting to note that the first three cumulants for T^/Tr — 10 at V g — 100r, 
200r, 250r, 300r and 400r have the same bias-voltagc-dcpendence as that for 
Yl/Tr = 0.1 at V g = 400r, 300r, 250r, 200r and 100r, respectively. In fact, 
besides the symmetry of the effective channel energy levels in Eqs. (|22|) and 
(123]) . Fermi distribution functions also satisfy the relation 

fL/R[4(i)}\ Vf - ed+u/2 ±AV = 1 _ /■H/i[ e i : (t)]|vi=e (t +l7/2 T AV ' ^ 
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fL/R[4(l)]\ Vg = ed+U/2 ±AV ~ 1 ~ [ £ 4 (t)] I V 9 =e d +t//2 T AV ' ( 25 ) 



which may ensure 



(m,2|p|2,m)| ri=ri rR=r2 = (to, 0| p |0, m)| r£=r2 rfl=ri , (26) 
± (m,l|p|l,™} ± | _ r = ± (m,l|p|l,m) ± | _ r r _ r - (27) 

1 L— L I: 1 K — 1 2 1 Z,— 1 2^ R — L 1 

where p (t) = ~J2 n (t). Based on these conditions one may prove the proba- 
bility distribution 

jP |2,m)| ri/=rii r R=r2 = P \0,m)\ rL=r2XR=ri , ( 2 8) 



= P + 
r -r r -r I 1 '" 1 ) 

1 t— l 111 H— 1 2 



(29) 

=r 2 ,r R =r, 



whose numerical results are shown in Fig. 2. It is because of the both symmetries 
of the effective channel energy levels and the probability distributions that the 
FCS of transport through the SMM for T L /T R = 10 at V g = e d + U/2 ± 
AV represents the same bias- voltage-dependence as that for T^/Tr = 0.1 at 
V g = €d + U/2 =p AV, which are shown in Fig. 1 for a given finite Coulomb 
interaction U — 100r. The transport characteristic should easily be observed 
experimentally by reversing the bias voltage between the left and right electrodes 
and tuning the applied gate voltage; but for the case of U — > oo, this feature 
will not be observed because the transition processes between doubly-occupied 
and the singly-occupied states are prohibited. 

It is generally thought that the Coulomb interaction gives rise to the sup- 



pression of shot noise 
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36] . For the present SMM system with finite Coulomb 
interaction, however, the super-Poissonian shot noise is observed in the sequen- 
tial tunneling regime. For a given U = 100r, as shown in Figs. 1(b) and (e), 
the super-Poissonian noise for Tl/Tr = 10 occurs at higher gate voltage such 
as V g = 300r and 400r, while that for Tl/Tr — 0.1 occurs at lower gate volt- 
age such as V g — 200r and 100r, which is consistent with the U — > oo caseQ. 
The presence of the super-Poissonian noise in the sequential tunneling regime 
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can be understood with the help of the dynamic competition between effective 

~nr r 



fast and slow channels 
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38, 39]. In order to give a qualitative 



explanation, we plot the six main molecular channel currents for V g — 400r as 
a function of bias voltage V& in Fig. 3(a). Here, the calculation of the molec- 
ular channel currents can be found in Refs. Q. When the bias voltage 

increases up to about 170r, the fast channel current I\2,2) ►11,5/2) begins to 

compete with the slow channel currents 1|2,— 2)->-|i,— 5/2) i -^|2,±2)->|i ±3/2)- an d 
^12 ±i)-y|i ±3/2)~ 1 but the competition is quickly destroyed due to the sum of slow 
channel currents being approach to the fast channel current, the corresponding 
noise [the short dashed line in Fig. 1(b)] is only enhanced but does not reach 
the super-Poissonian. With further increase of the bias voltage (about 250r), 
three sets of the fast-and-slow channel currents are formed, i.e., /|2,±2)->|i,±5/2) > 
^|2 ±2)->li ±3/2)~ ano - 1\2 ±i)-v|i ±3/2) ~ ) wn ich result in the second enhancement 
of shot noise. When the bias increases to about 500r, the fast channel cur- 
rents (J|2,2> — Hi.5/2), ^|2,2)->-|i,3/2>- and ^|2,i)-^|i,3/2)-) begin to decrease, while 

the corresponding slow channel currents (ija,— 2) >|i,-5/2) 5 I\2 -2)-s-|i -3/2) - and 

^|2 x> >-|x -3/2)~) begin to increase. The competition between the increase of 

the slow channel currents and the decrease of the fast channel currents finally 
leads to appearing of super-Poissonian Fano factor in the bias from about 500r 
to 625r [the short dashed line in Fig. 1(b)]. Furthermore, we also observe that 
the super-Poissonian distribution of the skewness F 3 in the sequential tunneling 
regime [Fig. 1(c)], which seems sensitive only to the competition between the 
fast channels of current decreasing and the slow channels of current increas- 
ing [Fig. 3(a)]. As for Tl/Tr = 0.1 at lower gate voltage, the mechanism of 
shot-noise enhancement originates from the same reason but the effective fast- 
and-slow transport channels consist of the transitions from the singly-occupied 
states to empty states, see Fig. 3(b). Furthermore, in the Coulomb blockade 
regime (see Fig. 1), the shot noise enhancement, as explained in Ref. jjyj], is 
due to the possible thermal occupation and subsequent sequential depletion of 
excited states that lead to small cascades of tunneling events interrupted by long 
Coulomb blockages. In fact, since in the Coulomb blockade regime the current 
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is exponentially suppressed and the electron transport is dominated by cotun- 
neling, when taking into account cotunneling the normalized second and third 
moments will deviate from the results obtained by only sequential tunneling 4C | . 



In contrast with the above case, we find that the finite Coulomb interaction 
plays a crucial role in determining whether the super-Poissonian noise occurs in 
situations with relatively small gate voltage for T l /Fr > 1 and relatively large 
gate voltage for Tl/Tr < 1. In order to study the effects of finite Coulomb 
interaction on the shot noise in transport through the SMM for the two above- 
mentioned cases, Fig. 4 shows the current and the shot noise as a function 
of the bias voltage for various Coulomb interaction energies. For the case of 
U — > oo, the super-Poissonian noise, as shown in Figs. 4(b) and (d), dose not 
appear in the sequential tunneling regime (also see Fig. 2 in Ref. [6[). Here, 
for the case of Tl/Tr = 0.1, U — > oo and V g — 300r, it should be noted 
that the super-Poissonian noise occurs in the Coulomb blockade regime, see the 
solid line in Fig. 4(d). Comparing with the U — > oo case, the finite Coulomb 
interaction may induce the super-Poissonian noise only when it is larger than a 
certain value. We take the case of V g = 100r and Tl/Fr = 10 for illustration. 
For small Coulomb interaction (for example, U = 100r) the super-Poissonian 
shot noise does not appear while for relatively larger Coulomb interaction (for 
example, U > 200r) the super-Poissonian shot noise is observed [see Fig. 4(b)]. 
The role of Coulomb interaction in enhancing shot noise can be understood 
with the help of the main channel currents shown in Fig. 5. The currents in 
Figs. 5(a) and (c) result from the electron transitions from singly occupied to 
empty states while Figs, (b) and (d) correspond to the transitions from doubly 
occupied to singly occupied states. For small Coulomb interaction U — 100r, 
when the bias changes from 200r to 500r the competition between the fast 

channel current (I\x t 5/2) ►|o,2)) an d the slow channel currents (/m -5/2) ►10,-2); 

/ij ±3/2)- y\o ±2) an d I\i ±3/2)" ►lo ±1)) * s an ac tive competition in which the 

increase (or decrease) of the fast channel current is always accompanied by the 
decrease (or increase) of the slow channel currents, which makes the shot noise 
increase from 0.84 to 0.93 [see Fig. 4(b)]. With the bias approaching to 600r the 
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currents I\ 2 ,2) — >\ifi/i) and I\ 2 2) *.|i >3 / 2 )± begin to compete with /|i,5/2>— >|o,2>, 

but the rapid increase of the amplitudes of I\2.2) >|i,5/2) an< ^ I\2 2) >\i3/2) ± 

lead to rapid decrease of -^1,5/2) >\o.2) and even destroy the competition, as a 

result the shot noise quickly drops to sub-Poissonian after increasing to 0.95 near 
Vb = 640r[see fig. 4(b)]. When Coulomb interaction increases to 200r, the bias 
voltage regime in which the active competition occurs is much larger than that 
for U = 100r, i.e., from 200r to 750r [see Fig. 5(c)], so that the shot noise has 
reached super-Poissonian before the effective competition between the fast and 

slow channels is destroyed by the rapid rise of I\2,2) Hi.5/2) an< ^ -^|2 2) Hi 3/2)* 

[see Figs. 5(c) and (d)]. 

On the other hand, although the shot noise of finite Coulomb interaction for 
r.L/r_R = 0.1 in Fig. 4(d) have the same bias- voltage-dependence as that for 
Fl/Fb. = 10 in Fig. 4(b), the corresponding transport processes are different 
due to occurring at different gate voltages, which may be found by comparing 
the main channel currents in Fig. 6(a) (for Tl/Yr = 0.1, U = 100r and 
V g = 400r) with those in Figs. 5(a) and (b) (for T L /T R = 10, U = 100r 
and V g = 100r). In particular, for the case of T L /r R = 0.1, U = 100r and 

V g = 400r there is a reverse current -T|o,2) Hi, 5/2) near Vb = 600r besides a 

rapidly increasing current >|o,2) [ see Fig- 6(a)] and the two currents 

destroy the competition between the fast and slow channels, resulting in rapid 
decrease of the shot noise. For the case of U — 200r and V g — 500r, as shown 
in Fig. 6(b), the bias voltage regime in which the active competition occurs is 
much larger than that for U = 100r and V g = A00T, i.e., from 200r to 750r (it 
is from 200r to 640r for U = 100r and V g = 400r), so that the shot noise has 
reached super-Poissonian before the reverse current and the rapidly increasing 
current occur. As a result the shot noise rapidly drops from the super-Poissonian 
to the sub-Poissonian. 

The theoretical investigations have predicted that the super-Poissonian noise 
may occur in the quantum dot systems 



38l . [39j, |40| and molecular junction 



13] 



coupled symmetrically to two electrodes. For the present SMM system, when it 
symmetrically couples to two leads the super-Poissonian noise in the sequential 
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tunneling regime is observed only under an appropriate finite Coulomb interac- 
tion that determined by the parameters of the SMM (see Fig. 7), which still not 
appears in the case of U = and U — > oo. The noise characteristics can still 
be attributed to the competition between the effective fast and slow transport 
channels. 

4. Conclusions 

We have studied the effect of finite Coulomb interaction on the FCS of elec- 
tron transport through a SMM weakly coupled to two metallic electrodes above 
the sequential tunneling threshold by means of the particle-number-resolved 
quantum master equation. The electron transport through the SMM with finite 
Coulomb interaction U has richer current correlation due to the transitions of 
doubly-occupied to singly-occupied states participating in transport. Compared 
with the U — > oo case, our analytical results show that the finite Coulomb in- 
teraction induces that the effective channel energy levels move symmetrically 
into the bias voltage window with increasing the bias voltage, as a result the 
FCS for V g = £d + U/2 ± AV has the same bias- voltage-dependence as that 
for V g = €d + U/2 =F AV (i.e., the FCS shows the symmetrical gate- voltage- 
dependence) when both the intensities of the SMM coupling to the left and 
right electrodes are interchanged, which should easily be observed experimen- 
tally by reversing the bias voltage. Moreover, we find that the effect of finite 
Coulomb interaction on the shot noise depends not only on the gate voltage, 
but also on the left-right asymmetry of SMM-clcctrode couplings. In the case 
of Tl/Th > 1, for an arbitrary given U (which contains U — and U — » oo) the 
super-Poissonian shot noise in the sequential tunneling regime may be observed 
at a relatively large gate voltage; whereas at a relatively small gate voltage the 
super-Poissonian shot noise does not occur for the cases of U = and U — > oo, 
and may appear only when U is the finite value which is related to the param- 
eters of the SMM. For the Tl/Tr < 1 case, the super-Poissonian shot noise 
is found at a relatively small gate voltage for an arbitrary given U; whereas 
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at a relatively large gate voltage it may occur only for a relatively large finite 
Coulomb interaction, which is contrary to the T^/Tr > 1 case. These charac- 
teristics of shot noise can be understood as a result of the active competition 
between the fast and slow channel currents. 
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Appendix 

In this appendix, we take (l, — 1 1 S (x, t) | 1, — § ) as an example to illustrate 
the procedure for obtaining the specific form of C x . Based on the S (%, i) def- 
inition and Eq. (|11[) . after a careful calculation, the equation of motion for 



<1, 



S {X,t) |l, -§) is given by 




€ 



(2,-2|£|2,-2) 



+T R n 



(- 

R 



2,-2) -%-$) 



e lx (2,-2|5|2,-2) 




- e |0,-2> 



(0,-2|5|0,-2) 



+T R n 



(+ 

R 



e |l,-|) " £ |0,-2> 



e' tx (0, -21510, -2). 



(A.1) 
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Similarly, we can carry out the equations of motion for other matrix elements, 
which can be rewritten as a compact matrix form 



S = C X S. 



Here, the column matrix S has the following form 

5' 



(A.2) 



S 



1. 



s 



1-- 



<2,-2|S|2,-2), 



(0,-2|£|Q,-2) +(1, 



S 



1, 



,5' 



L — 



1,-: 



(0,-l|5|0,-l),+ (l, 



S 



1, 



,5' 



1,— 



1,-: 



1,-: 



,(2,-l|S|2,-l) 



1,-7 



1,-: 



1,-: 



,(2,0|S|2,0), 



(0,0|5|0,0) 


s 




s 






s 




(4 


s 


1^) ,{2,11512,1}, 


(0,1|S|0,1) + 


s 




s 


'■!>"• 



' 2 



S 



1, 



1, ■ 



'2 



<2,2|£|2,2) 



(0,21510,2), (1, 



From Eq. (|A.1[) . the first row of C x is given by 

[^x]lx28 ~ 



[c x \ !2 ,[c x l 3 , 0,0, 0,0, 0,0, 

0,0,0,0,0,0,0,0, 

0,0,0,0,0,0,0,0, 

0,0,0] 



(A.3) 



(A.4) 
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with 

[^xll,! = ~ r i"i +) ( £ |2,-2) - £ |l,-f }) ~ T R n R ) ( £ |2,-2) - £ |l,-f }) 

-Tlu^ (e|i,_|) -e|0,-2>) - Trti^ (c|i,_j) - £|o,-2>) , 

[Atli,2 = r i n i _) ( e |2,-2> - e |i,-|)) + (e| 2 ,-2> - e |i,-|)) e ' X : 

[ £ x]i, 3 = T Ln ( L ] (e|i,_|) - e|0,-2>) + !V4 +) (e|i,_|) - C|o,-2)) e~* x - 
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Figure 1: (Colour online) The first three cumulants of zero-frequency current fluctuation 
versus bias voltage for different gate voltages, (a), (b) and (c) for Tl/Tr = 10, (d), (e) and 
(f) for T L /T R = 0.1. The molecular parameters: S = 2, e d = 200r, J = 100I\ K 2 = 40r, 
B = 80r and k B T = WF. 
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Figure 2: (Colour online) The probability distribution of molecular eigenstates versus bias 
voltage for different gate voltages with U = 100r. (a), (b) and (c) V g = 4001" and r L /T R = 10, 
(d), (e) and (f) V g = lOOT 1 and Tl/Tr = 0.1. The molecular parameters are the same as in 
Fig. 1. 
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Figure 3: (Colour online) The channel currents versus bias voltage for various gate voltages 
with U = 100r. (a) V g = 400r and T L /T R = 10, (b) V g = 100V and V L /r R = 0.1. The 
\i) — > \j) denotes the molecular channel current The molecular parameters are the 

same as in Fig. 1. 
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Figure 4: (Colour online) The average current and shot noise versus bias voltage for the 
different Coulomb interaction energies, (a) and (b) for T^/Tr = 10, (c) and (d) for Tl/Tr = 
0.1. The molecular parameters are the same as in Fig. 1. 
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Figure 5: (Colour online) The channel currents versus bias voltage for different Coulomb 
interaction energies with T L /T R = 10. (a) and (b) U = WOT, V g = 100F, (c) and (d) 
U = 200r, V g = lOOr. The channel current expression and molecular parameters are the 
same as Fig. 3. 
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Figure 6: (Colour online) The channel currents versus bias voltage for different Coulomb 
interaction energies with Y L /T R = O.f. (a) for U = 100r, V g = 400I\ (b) U = 200I\ 
V g = 500r. The channel current expression and molecular parameters are the same as Fig. 3. 
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Figure 7: (Colour online) The average current and shot noise versus bias voltage for different 
Coulomb interaction energies with = Tr = T. The molecular parameters are the same as 
in Fig. 1. 



